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ABSTRACT 

We report on the X-ray and multiwavelength properties of 1 1 radio-quiet quasars with weak or no emission 
lines identified by the Sloan Digital Sky Survey (SDSS) with redshift z = 0.4-2.5. Our sample was selected 
from the Plotkin et al. catalog of radio-quiet, weak-featured AGNs. The distribution of relative X-ray brightness 
for our low-redshift weak-line quasar (WLQ) candidates is significantly different from that of typical radio- 
quiet quasars, having an excess of X-ray weak sources, but it is consistent with that of high-redshift WLQs. 
Over half of the low-redshift WLQ candidates are X-ray weak by a factor of > 5, compared to a typical SDSS 
quasar with similar UV/optical luminosity. These X-ray weak sources generally show similar UV emission- 
line properties to those of the X-ray weak quasar PHL 1811 (weak and blueshifted high-ionization lines, weak 
semi-forbidden lines, and strong UV Fe emission); they may belong to the notable class of PHL 1811 analogs. 
The average X-ray spectrum of these sources is somewhat harder than that of typical radio-quiet quasars. 
Several other low-redshift WLQ candidates have normal ratios of X-ray-to-optical/UV flux, and their average 
X-ray spectral properties are also similar to those of typical radio-quiet quasars. The X-ray weak and X-ray 
normal WLQ candidates may belong to the same subset of quasars having high-ionization "shielding gas" 
covering most of the wind-dominated broad emission-line region, but be viewed at different inclinations. The 
mid-infrared-to-X-ray spectral energy distributions (SEDs) of these sources are generally consistent with those 
of typical SDSS quasars, showing that they are not likely to be BL Lac objects with relativistically boosted 
continua and diluted emission lines. The mid-infrared-to-UV SEDs of most radio-quiet weak-featured AGNs 
without sensitive X-ray coverage (34 objects) are also consistent with those of typical SDSS quasars. However, 
one source in our X-ray observed sample is remarkably strong in X-rays, indicating that a small fraction of 
low-redshift WLQ candidates may actually be BL Lacs residing in the radio-faint tail of the BL Lac population. 
We also investigate universal selection criteria for WLQs over a wide range of redshift, finding that it is not 
possible to select WLQ candidates in a fully consistent way using different prominent emission lines (e.g., 
Lya, C IV, Mg II, and Hj3) as a function of redshift. 

Subject headings: galaxies: active — galaxies: nuclei — quasars: emission lines — X-rays: galaxies — BL 
Lacertae Objects: general 



1. INTRODUCTION 

Strong and broad line emission is a common feature of 
quasar spectra in the optical and UV bands. However, since 
multi-color quasar selection at high redshift in the Sloan Dig- 
ital Sky Survey (SDSS; York et al. 2000) is mostly based 
upon the presence of the Lya forest and Lyman break (e.g., 
Richards et al. 2002), the SDSS can also effectively select 
high-redshift quasars with weak or no emission lines. About 
90 such weak-line quasars (WLQs) at high redshift have 
been found with Lya H- N V rest-frame equivalent widths of 
REW < 15 A (e.g.. Fan et al. 1999, 2006; Anderson et al. 
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2001; Coflinge et al. 2005; Diamond-Stanic et al. 2009, here- 
after DS09). Some of these objects show a hint of weak Lya 
emission but no other lines; others are completely bereft of 
detectable emission lines even in high-quality spectra. High- 
redshift SDSS quasars show an approximately log-normal 
distribution of Lya H- N V REW with a mean of ss 62 A 
(DS09). The WLQs constitute >3c7 negative deviations from 
the mean, and there is no corresponding population with 
>3(7 positive deviations. The majority of these high-redshift 
WLQs are radio quiet (o!ro > —0.21; (Xto is the slope of a nom- 
inal power law between 5 GHz and 2500 A in the rest frame; 
see 21 for a full definition). 

WLQs have mainly been studied at high redshifts due to the 
fact that the Lya forest enters into the SDSS spectroscopic 
coverage for quasars at z > 2.2. However, there is no appar- 
ent reason to believe that these objects should not also exist 
at lower redshifts. Indeed, a few apparent analogs of WLQs 
at lower redshifts have been found serendipitously over the 
past w 15 years; e.g., PG 1407H-265 (McDowell et al. 1995; 
z = 0.94), 2QZJ2154-3056 (Londish et al. 2004; z = 0.49), 
and PHL 1811 (Leighly et al. 2007ab; z = 0.19). As a byprod- 
uct of a systematic survey for optically selected BL Lacertae 
objects (hereafter BL Lacs) in SDSS Data Release 7 (DR7; 
Abazajian et al. 2009), Plotkin et al. (2010a) discovered about 
60 additional radio-quiet WLQ candidates at z < 2 . 2 for which 
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all emission features have REW < 5 A. These objects are 
perhaps the first low-redshift SDSS counterparts of the pre- 
viously identified high-redshift SDSS WLQs. Following the 
nomenclature that has been established by previous work on 
WLQs (e.g., Shemmer et al. 2009), we define "high redshift" 
as z > 2.2 and "low-redshift" as z ^ 2.2, because WLQs are 
selected with different approaches for these redshift ranges 
(see above). Although WLQs are rare, their exceptional char- 
acteristics constitute a challenge to our overall understanding 
of quasar geometry and physics, especially the quasar broad 
emission-line region (BELR). Analogously, physical insights 
have been gained by investigating other minority populations 
with exceptional emission-line or absorption-line properties, 
such as Narrow-Line Seyfert 1 (NLSl) galaxies and Broad 
Absorption Line (BAL) quasars. Therefore, extensive studies 
of the multi-band properties of WLQs should have scientific 
value. 

There are several candidate explanations for the physical 
nature of WLQs. Their UV emission lines may be weak due to 
an "anemic" BELR with a significant deficit of line-emitting 
gas (e.g., Shemmer et al. 2010). It has also been speculated 
that WLQs may represent an early stage of quasar evolution in 
which an accretion disk has formed and emits a typical contin- 
uum, but BELR formation is still in progress (e.g., Hryniewicz 
et al. 2010; Liu & Zhang 201 1). 

The weak UV emission lines may also be a consequence of 
a spectral energy distribution (SED) which lacks high-energy 
ionizing photons. This soft SED may be a result of unusual 
accretion rate. For example, an extremely high accretion rate 
might produce a UV-peaked SED (e.g., Leighly et al. 2007). 
In this scenario, high-ionization lines, like C IV, should be 
suppressed relative to low-ionization lines like Hj3 . However, 
Shemmer et al. (2010) estimated the normalized accretion 
rates, L/Leaa, of two high-redshift WLQs via near-infrared 
spectroscopy and found their accretion rates were within the 
range for typical quasars with similar luminosities and red- 
shifts. Alternatively, a combination of low accretion rate and 
large black hole mass may lead to a relatively cold accretion 
disk that emits few ionizing photons. Laor & Davis (2011) 
predicted a steeply falling SED at A < 1000 A for quasars 
with cold accretion disks, and such an SED was observed in 
the WLQ SDSS J0945H-1009 by Hryniewicz et al. (2010). 

High-energy ionizing photons (including X-rays) may be 
heavily absorbed before they reach the BELR. Wu et al. 
(2011) studied a population of X-ray weak quasars with un- 
usual UV emission-line properties like those of PHL 1811 
(weak and highly blueshifted high-ionization lines, weak 
semi-forbidden lines, and strong UV Fe emission). All of 
their radio-quiet PHL 1811 analogs were found to be X-ray 
weak by a factor of w 13 on average. These objects also 
show a harder average X-ray spectrum than those for typ- 
ical quasars which suggests the presence of X-ray absorp- 
tion. PHL 1811 analogs appear observationally to be a sig- 
nificant subset (w 30%) of WLQs. The existence of a class of 
quasars with high-ionization "shielding gas" covering most 
of the BELR, but little more than the BELR, could potentially 
unify the PHL 1811 analogs and WLQs via orientation effects 
(see §4.6 of Wu et al. 2011). The shielding gas would absorb 
high-energy ionizing photons before they reach the BELR, re- 
sulting in weak high-ionization emission lines. When such a 
quasar is observed through the BELR and the shielding gas, 
a PHL 1811 analog would be seen; when it is observed along 
other directions, an X-ray normal WLQ would be observed. 



Another possibility is that instead of being intrinsically 
weak, the UV emission lines of WLQs could in principle be 
diluted by a relativistically boosted UV/optical continuum as 
for BL Lac objects. However, this scenario is not likely for 
most WLQs. Shemmer et al. (2009) found that the X-ray 
properties of high-redshift WLQs are inconsistent with those 
of BL Lac objects. Furthermore, there is no evidence of strong 
optical variability or polarization for these WLQs (see DS09; 
Meusinger et al. 2011). The UV-to-infrared SEDs of high- 
redshift WLQs are also similar to those of typical quasars, 
while the SEDs of BL Lac objects are much different (DS09; 
Lane et al. 2011). Nevertheless, it is possible that the pop- 
ulation of BL Lac objects has a small radio-quiet tail (e.g., 
Plotkin et al. 2010b) and that a small fraction (< 5%; see 
Lane et al. 2011) of the general WLQ population may be 
BL Lac objects. 

Most previous studies of WLQs were based on high- 
redshift objects. To investigate the nature of the overall 
WLQ population, we obtained new X-ray observations of 
low-redshift WLQs selected mainly from the catalog of radio- 
quiet BL Lac candidates in Plotkin et al. (2010a). We also 
utilized sensitive archival X-ray coverage of the sources in 
their catalog. Our closely related science goals are the fol- 
lowing: (1) enable comparison of the broad-band SEDs of 
low-redshift WLQs to those of high-redshift WLQs, typical 
radio-quiet quasars, and BL Lac objects; (2) provide basic 
constraints upon X-ray spectral properties via band-ratio anal- 
ysis and joint spectral fitting; (3) clarify if there is broad-band 
SED diversity among low-redshift WLQs; and (4) allow re- 
liable planning of future long, spectroscopic X-ray observa- 
tions. 

In f|2] we describe the selection of our sample of 
low-redshift, radio-quiet WLQ candidates. In ij3] we detail 
their UV/optical observations and the measurement of their 
rest-frame UV spectral properties. In 21 we describe the rel- 
evant X-ray data analyses. Overall results and associated dis- 
cussion are presented in f|5] Throughout this paper, we adopt 
a cosmology with Hq = 70.5 km s"' Mpc"', Q.^ = 0.274, 
and iiA = 0.726 (e.g., Komatsu et al. 2009). 

2. SELECTION OF THE LOW-REDSHIFT WLQ 
CANDIDATES 

We obtained Chandra snapshot observations (3.0-4.1 ks) 
of six low-redshift (z = 0.40-1.67) WLQ candidates. Five 
of the six targets were identified by Plotkin et al. (2010a) 
as radio-quiet, weak-featured SDSS quasars with all emis- 
sion features having REW < 5 A. An additional source, 
SDSS J0945H-1009, was similarly identified as a weak- 
featured quasar by Hryniewicz et al. (2010). All the objects 
are sufficiently bright in the optical band (m, < 18) for short 
Chandra observations to provide tight constraints on their X- 
ray-to-optical SEDs. 

We further utilized the weak-featured quasar catalogs in 
Plotkin et al. (2010a) to search for low-redshift, radio-quiet 
sources having sensitive archival X-ray coverage. To en- 
sure our sample has the high X-ray detection fraction nec- 
essary to provide physically meaningful constraints, we only 
selected sources covered by Chandra or XMM-Newton obser- 
vations An additional five sources were thereby added into 
our sample. Three of them (J1013H-4927, Jl 139-0201, and 

' We also checked for pointed ROSATFSPC observations witli an exposure 
time greater than 5 ks and an off-axis angle less than 19' (within the inner ring 
of the PSPC detector). However, none of the radio-quiet, low-redshift sources 



WEAK EMISSION-LINE QUASARS AT LOW REDSHIFT 



3 




» Chondro Cycle 12 torgets 

k Archival sources from Plotkin et ol. (2010Qb) 

I RQ WLOs from Shemmer et ol, (2006, 2009) 



Fig. 1. — SDSS absolute ('-band magnitude, Mj, plotted versus redshift, z- 
The red filled circles and triangles show our sample of low-redshift WLQ 
candidates; the blue filled squares show high-redshift WLQs from Shemmer 
et al. (2006, 2009); the grey dots represent the 105,783 objects in the SDSS 
DR7 quasar catalog (Schneider et al. 2010). 

J 1604+4326) appear in the radio-quiet, weak-featured quasar 
catalog (Table 6 in Plotkin et al. 2010a). Jl 139-0201 was tar- 
geted by Chandra as an optically selected BL Lac candidate in 
Cycle 5, while J1013+4927 and J1604+4326 were serendip- 
itously covered by Chandra or XMM-Newton observations. 
The other two objects (J21 15+0001 and J2324+1443) were 
initially identified as weak-featured quasars by Collinge et al. 
(2005). They were also listed in the catalog of Plotkin et al. 
(2010a). These two sources did not have constraints on their 
radio fluxes in Collinge et al. (2005) or Plotkin et al. (2010a) 
but were later confirmed as radio-quiet sources by the VLA 
observations of Plotkin et al. (2010b). They were targeted by 
Chandra as radio-quiet BL Lac candidates in Cycle 10; their 
observations were briefly reported in Plotkin et al. (2010b). 
Table [Tlpresents the X-ray observation log for our sample. 

Our sample includes 1 1 WLQs in total. All of the sources 
in our sample have redshifts of z < 2.2, except J21 15+0001 
which has a slightly higher redshift of z = 2.4995 (see 
for redshift measurements). For comparison, all the radio- 
quiet WLQs studied in X-rays by Shemmer et al. (2006, 2009) 
have z > 2.7 (see Fig. Fig. [2] shows the SDSS spectra 
of the sources in our sample. The spectra show no evidence 
for dust reddening or intrinsic broad absorption lines (BALs); 
i.e., there is no indication that their UV/optical continua or 
BELRs are obscured. We will compare the multiwavelength 
properties of our sample to those of the high-redshift WLQs 
in Shemmer et al. (2006, 2009) in §|5] 

3. UV/OPTICAL OBSERVATIONS 

3.1. UV Emission-Line Measurements 

The redshift values (see Table[T]i for our low-redshift WLQs 
are generally those from Hewett & Wild (2010) which are 
the best available measurements for large SDSS quasar sam- 
ples. There are three sources lacking Hewett & Wild mea- 
surements. For two quasars (Jl 109+3736 and Jl 139— 0201), 
the redshift values are taken from the catalog of Plotkin et al. 
(2010a). The redshift of the other source (J21 15+0001; z = 
2.4995 ± 0.0052) is measured based on a Ly a + C IV absorp- 

in the catalogs of Plotkin et al. (2010a) is covered by ROSAT observations 
meeting these criteria. 



tion system[3 

To obtain accurate measurements of the weak emission 
lines, we manually measured rest-frame emission-line prop- 
erties for C IV, Si IV, the A 1900 complex^ and Fe III UV48 
(see Table|2]i following the method in §2.2 of Wu et al. (201 1), 
which is summarized below. We first smoothed the SDSS 
spectra with a 5-pixel sliding-box filter, and manually inter- 
polated over strong narrow absorption regions. We then fit- 
ted a power-law local continuum for each line between their 
lower and upper wavelength limits Aio and Ahi (see Table 2 
of Vanden Berk et al. 2001). After subtracting the local 
continuum, we measured the REW value for each line. The 
C IV blueshifts were calculated between the lab wavelength 
in the quasar rest frame (1549.06 A, see Table 2 of Vanden 
Berk et al. 2001) and the observed mode of all pixels with 
heights greater than 50% of the peak height, where mode = 
3 X median— 2 X mean. For comparison, we also include in 
Table |2] the corresponding measurements of the spectrum of 
PHL 1811 (Leighly et al. 2007b) and of the composite spec- 
trum of typical SDSS quasars in Vanden Berk et al. (2001). 
The spectral measurements of PHL 1811 are included here be- 
cause some of our low-redshift WLQ candidates show similar 
unus ual U V/optical spectral properties to those of PHL 1811 
(see i j5.2l i. The Mg II measurements from Shen et al. (2011) 
are also listed in Table |2] These measurements are reliable 
because the Fe II component, which could affect the Mg II 
strength measurement, was well modeled. These REW(Mg ll) 
values somewhat exceed the selection criterion of REW < 5 A 
for BL Lac candidates in Plotkin et al. (2010a). This dis- 
crepancy mainly originates from differences in measurement 
methods. For Plotkin et al. (2010a), it was impractical to 
define reference wavelengths to model the continuum in a 
uniform way for the entire large sample since many objects 
lack redshift measurements. The REW values in Plotkin et al. 
(2010a) were measured manually after defining the contin- 
uum by eye for most sources. While this method generally 
performed well for BL Lac objects, it did not properly model 
blended Fe emission for unbeamed objects. 

Only two sources (J0945+1009 and J 16 12+5 118) have 
high-quality C IV coverage in their SDSS spectra so that 
we are able to measure their C IV REW and blueshift val- 
ues. Both sources have weak and highly blueshifted C IV 
lines. Jl 139— 0021 has no clearly detectable C IV line in its 
SDSS spectrum; we could only obtain an upper limit upon its 
REW. Therefore, we obtained follow-up UV spectroscopy for 
this source with the Low-Resolution Spectrograph (LRS; Hill 
et al. 1998) on the Hobby-Eberly Telescope (HET; Ramsey 
et al. 1998). The UV emission-line measurements based on 
the HET spectroscopy are also listed in Table|2] Jl 139— 0021 
has a weak and strongly blueshifted C IV line in its HET spec- 
trum. 

All of the sources having C III] coverage show weaker C III] 
semi-forbidden lines than those of typical quasars. The Fe III 
UV48 strength of our low-redshift WLQ candidates is gener- 
ally similar to those of typical quasars. The SDSS spectrum of 

Plotkin et al. (2010b) did not report the redshift for this source. In this 
work we adopt the redshift of the Lya + C IV narrow absorption system as 
the systemic redshift. Nestor et al. (2008) fit a Gaussian distribution centered 
at V = km s^' with a = 450 km s^' to the distribution of narrow C IV 
systems around quasar systemic redshifts. We measured the redshift using 
that Gaussian dispersion as the redshift uncertainty to obtain j = 2.4995 ± 
0.0052. 

" Mainly C III] A 1909, but also including other features; see Note (b) of 
Tabled 
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Fig. 2. — SDSS spectra for the 1 1 sources in our sample of low-redshift WLQ candidates, ordered by AcHqx (see i|4]for definition). The AOox values and their 
error bars (if the source is detected in X-rays) are shown for each source. The name of each source is labeled in the format of 'ihhmm+ddmni . The i -coordinates 
are the flux density (Fi) in arbitrary linear units. The tick marks on the y-axis show the zero flux-density level for each normalized spectrum. The spectra have 
been smoothed using a 5-pixel sliding-box filter The spectrum of the radio-quiet BL Lac candidate Jl 109-H3736 is shown separately in the lower right panel for 
convenience of presentation (since its redshift is much lower than those of the other sources in our sample). Emission fines, including C IV A 1549, C III] A 1909, 
and Mg II A2799, are labeled in the left and upper right panels. The Mg II A2799, H/3 A4862, [O III] A5007, and Ha A6564 fines are labeled in the lower right 
panel; the Ca II H/K break is also marked by the dotted lines. All the quoted values here are vacuum wavelengths. The spectral resolution is, 2000. Also 
included are the composite spectrum of SDSS quasars by Vanden Berk et al. (2001) and the mean spectrum of the high-redshift WLQs of Shemmer et al. (2009). 
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Jl 109+3736 does not have coverage of these rest-frame UV 
emission lines because of its much lower redshift, while the 
signal-to-noise ratio of the SDSS spectrum of J21 15+0001 is 
too low to make reliable emission-line measurements. 

3.2. Comparing the Emission-Line Strengths of 
Low-Redshift WLQ Candidates and Typical SDSS 
Quasars 

After measuring the strengths of the UV emission lines of 
our low-redshift WLQ candidates (see Table|2]i, we further in- 
vestigated the REW distributions of prominent emission lines 
that are covered by the SDSS spectra of low-redshift quasars 
(such as C IV, Mg II, and H/3). This allowed assessment of the 
emission-line weakness of our low-redshift, radio-quiet WLQ 
candidates compared to typical SDSS quasars. Furthermore, 
our method of selecting low-redshift WLQ candidates is dif- 
ferent from that for high-redshift WLQs because low-redshift 
quasars do not have coverage of Lya + N V emission in their 
SDSS spectra. The REW distributions of emission lines like 
C IV, Mg II, or H/3 may give insights on universal selection 
criteria for WLQ candidates at different redshifts. 

We utilized the REW measurements of C IV, Mg II, and 
Hj3 from the catalog of Shen et al. (2011) for SDSS DR7 
quasars (see their Table 1). For the C IV and Mg II lines, we 
adopted the REW values for the entire lines, while for Hj8 
we added the REW values of the broad and narrow compo- 
nents to obtain the REW of the entire line. Shen et al. (201 1) 
reported the REW distributions of these lines for all SDSS 
DR7 quasars having applicable REW measurements (see their 
Figs. 12-14). In this work, we selected unbiased samples of 
SDSS quasars with high-quality optical/UV spectra to study 
the REW distributions of these emission lines by imposing the 
following criteria: 

1. We only use the DR7 quasars selected with the final 
algorithm given by Richards et al. (2002) to maintain 
consistency with DS09. 

2. BAL quasars that were cataloged in Gibson et al. 
(2009) and Shen et al. (201 1) were removed. 

3. We restricted the redshift ranges of the objects for each 
line (C IV: z = 1.55^.67; Mg II: z = 0.42-2.15; Hj3: 
z < 0.84) to ensure that the SDSS spectra of these ob- 
jects cover the whole region of each emission line de- 
fined by Aio and Ay listed in Table 2 of Vanden Berk 
etal. (2001). 

4. To select objects with high-quality SDSS spectra, we 
calculated the signal-to-noise ratio (S/N) of continuum 
regions close to each emission line, SNnoo for C IV, 
5A^30oo for Mg II, and 5A^5i5o for Hj3 (see Tabled, fol- 
lowing the method of Gibson et al. (2009). SNnoo, 
SNiQoo, and SNsiso are calculated as the median of the 
ratio between the flux and the error (obtained from the 
SDSS pipeline) for all the spectral bins in the rest-frame 
1650-1750 A, 2950-3050 A, and 5100-5200 A re- 
gions, respectively. These wavelength regions are free 
of strong emission and/or absorption features, but are 
still close to the above emission lines in our study. We 
require each of the above continuum S/N to be greater 
than 7. The S/N values of the emission lines them- 
selves were not utilized because that would introduce 
bias against objects with weak emission Unes. 



5. We eliminated the objects that have large fractions of 
bad pixels in their SDSS spectra in the wavelength 
ranges of these emission lines (these lead to unreliable 
measurements). We imposed the following cuts on the 
numbers of pixels that were included in the fitting for 
each line given in the catalog of Shen et al. (2011): 
LINE_NPIX_CIV > 250 for C IV, LINE_NPIX_MGII 
> 300 for Mg II, and LINE_NPIX_HB > 150 for Hj3. 

Applying the above cuts on redshift, continuum S/N, and the 
numbers of pixels in the fitting (criteria 3-5 above) removed 
35%, 27%, and 37% of the objects in the REW distribution 
investigations for C IV, Mg II, and Hj3, respectively. All the 
above restrictions and quality cuts are necessary because un- 
reliable line measurements could significantly affect the REW 
distributions particularly in the tails with low or high REW 
values. 

Fig.[3]shows the REW distributions of the C IV, Mg II, and 
Hj3 lines for our selected samples of SDSS quasars. We also 
include the histogram for Lya + N V from DS09 for com- 
parison. DS09 also showed a histogram for C IV REWs of 
SDSS DR5 quasars, which has a similar profile as the C IV 
REW histogram in Fig. |3] We fit the REW histogram of each 
line with a lognormal distribution (see the blue solid lines in 
Fig.|3l also see the dotted and dashed lines for the 2(7 and 3(7 
ranges of each lognormal model). The REW measurements 
for the quasars in our low-redshift WLQ sample (see 
are also shown in Fig. [3] All the C IV REW values of our 
sources are far below the negative 3(7 deviation of the log- 
normal distribution. Most of the REW(Mg ll) values are also 
below the negative 3(7 deviation of the lognormal distribution; 
the largest REW(Mg ll) value for our sample is close to the 
negative 2(7 deviation (see the top right panel of Fig. [3]). 

The REW distribution of Lya + N V in DS09 shows a 
prominent tail toward low REW values; this tail is the basis 
on which the high-redshift WLQs are defined. The C IV REW 
distribution shows similar behavior in that a more prominent 
skew tail toward low REW values exists, while there is no 
corresponding tail toward high REW values. However, the 
histogram of the Mg II REWs appears more symmetric, with 
only small tails toward both the low end and high end of the 
REW distribution. The Hj3 REW distribution is similar to 
that for Lya + N V; a prominent tail toward low REW val- 
ues exists, while the tail toward high REW values is much 
less significant. We randomly chose sets of sources in the 
tails with low or high REW values in the histograms for C IV, 
Mg II, and Hj3 and then visually examined their SDSS spectra 
and the quality assessment plotQfor their individual spectral 
fits in Shen et al. (2011). These sources have good spectral- 
fit quality; their REW measurements should be reliable. It 
is worth noting that the quasars in the REW histograms for 
various lines have different ranges of redshift and luminosity, 
which may affect their REW distributions (e.g., via the Bald- 
win effect; Baldwin 1977). We test this hypothesis by com- 
paring the C IV and Mg II REW distributions for a set of SDSS 
quasars with 1.55 < z < 2.15 and 45 < logvLgooo < 46 
(VL3000 is the luminosity at rest-frame 3000 A in erg s^^ ob- 
tained from Shen et al. 2011). Fig. |4] shows that both the 
C IV and Mg II REW distributions show similar behavior to 
the distributions presented in Fig [3j the C IV histogram has a 
significant tail at the low REW end, while the Mg II histogram 
is symmetric with weak tails. While the physical reason for 

See https://www.cfa.harvard.edu/ yshen/BH_mass/dr7.htm. 
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Fig. 3. — Distributions of REWs of Lya + N V, C IV, Mg II, and H/3 for SDSS DR7 quasar samples described in i|3.2l In each panel, the blue solid line shows 
the best-fit lognormal distribution. The dotted and dashed lines show the 2f7 and 3f7 ranges, respectively. The 3f7 range is also noted in the upper right corner of 
each panel. The red filled squares with arbitrary y-coordinates are the REW values for our current sample. The redshift range and the number of sources for each 
line REW distribution are also noted. 



the different emission-line REW distributions is unclear, the 
Mg II line has low optical depth, while the C IV line has much 
higher optical depth (e.g., see Eracleous et al. 2009). In the 
context of a disk-wind model for the BELR (e.g., Murray & 
Chiang 1997; Proga et al. 2000), the C IV emission is consid- 
ered to be mainly from the AGN wind, while the Mg II line 
mostly originates from the accretion disk (e.g., Leighly 2004; 
Richards etal. 2011). 

The 3(7 tail toward low REW values of the C IV histogram 
is defined by REW(C IV) < 10 A, which is consistent with 
that derived in DS09. Using this criterion one could perhaps 
hope to extend the redshift range of WLQ selection from z > 
3 down to z > 1 .5. However, this definition of WLQs has sig- 
nificant inconsistency with that in DS09 based on REW(Lya 
-I- N V) < 15.4 A. Fig.|5jfl) shows the relation between the 
REW values of Lya H- N V and C IV, and the best-fit power- 
law model found using the IDL routine LINMIX.ERrQ The 

The LINMIX_ERR procedure is a Bayesian approach to linear regression 
which usually has good performance when there is significant intrinsic scatter 
and con'elated eiTor bars (see Kelly 2007 for more details). To identify the 
best-fit power-law model, we first fit the correlation by assigning REW(Lyo: 
+ N V) as the "independent" variable and REW(C IV) as the "dependent" 
variable, and then exchange these two variables to obtain another fitting cor- 
relation. We finally calculated the bisector of these two power-law models as 



shaded area shows the 90% confidence uncertainty range ob- 
tained via a nonparametric bootstrap method (Efron 1979). 
Although the REW values of these two lines are positively 
correlated, large scatter exists. Over half of the sources with 
REW(Lya h- N v) < 15.4 A have REW(C iv) > 10 A, and 
vice versa (see Fig.|5]a!)- The large scatter is likely to be intrin- 
sic, since the REW measurement errors (see Fig.|5]for typical 
error bars) are much smaller than the scatter DS09 also sug- 
gested the inconsistency between sources with weak Lya -i- 
N V and those with weak C IV emission. A total of 39 of 
the 74 WLQs cataloged in DS09 have REW(C IV) > 10 A 
based on the measurements in Shen et al. (201 1)13 As stated 
in DS09, their measurements may underestimate REW(Lya 
-I- N V) when strong intervening absorption exists; this is per- 
haps one reason for the inconsistency. The correlation be- 
tween the REWs of C IV and Mg II also has significant scat- 
ter, and so does the correlation between the REWs of Mg II 
and Hj3 (see Figs.[2> and|5]:). The red shaded histograms in 
Fig.|4]show the REW(C IV) distribution for sources with weak 

the best-fit model (e.g., Isobe et al. 1990). We used the same bisector method 
for the correlations between the REW values of other emission fines. 

It is worth noting that many WLQs in the DS09 sample have large un- 
certainties on their C IV measurements. Only 13 of the 74 WLQs definitely 
have REW(C IV) > 10 A at 3(T significance. 
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Fig. 4. — Distributions of REWs for C IV (top panel) and Mg II (bot- 
tom panel) for the same set of SDSS quasars with 1.55 < z < 2.15 and 
45 < logLsooo < 46. The red shaded histogram in the top panel shows the 
REW(C IV) distribution for sources with REW(Mg II) below the 3fT negative 
deviation for Mg II. The red shaded histogram in the bottom panel shows the 
REW(Mg II) distribution for sources with REW(C IV) below the 3a negative 
deviation for C IV. The bottom panel shows that many objects with weak C IV 
emission do not have weak Mg II. Other lines follow the same definitions as 
those in Fig.|3] 

Mg II and the REW(Mg ll) distribution for sources with weak 
C IV. While the sources with weak Mg II also tend to have 
weak C IV (below the negative 2(7 deviation), many objects 
with weak C IV have fairly strong Mg II (see the bottom panel 
of Fig.©. 

Given the results above, it is therefore difficult to find 
consistent criteria for WLQs at different redshifts using dif- 
ferent emission lines even though the REW distributions of 
C IV and Hj3 show similar behavior to that of Lya -i- N V. 
Since there is no single line that is covered by SDSS spec- 
troscopy for quasars at all redshifts between zero and six, it 
appears that there is not a straightforward, direct way to define 
universal selection criteria for WLQs at all redshifts solely 
based on SDSS spectroscopy. Therefore, we only choose 
our low-redshift WLQ candidates mainly from the catalog 
of Plotkin et al. (2010a), which has a strict criterion on 
all emission-line strengths (REW < 5 A). Future UV spec- 
troscopy which covers the Lya -i- N V and/or the C IV regions 
for low-redshift WLQs will provide insights toward a univer- 
sal definition for WLQs. 

4. X-RAY DATA ANALYSIS 

The six new low-redshift WLQ candidates targeted in 
Chandra Cycle 12 were observed with the S3 CCD of 
the Advanced CCD Imaging Spectrometer (ACIS; Garmire 
et al. 2003). The reduction of the Chandra data was per- 
formed using standard CIAO v4.3 routines. X-ray images 
were produced for the observed-frame soft (0.5-2.0 keV), 
hard (2.0-8.0 keV), and full (0.5-8.0 keV) bands using ASCA 
grade 0, 2, 3, 4, and 6 events. The WAVDETECT algorithm 
(Freeman et al. 2002) was run on the images using a de- 
tection threshold of 10^^ and wavelet scales of 1, \/2, 2, 
2^2, and 4 pixels. AU targets, except J0945+1009, were 



detected by Chandra within 0.8" of the optical coordinates. 
The X-ray images of J0945+1009 were visually examined, 
and no hint of a detection was found. Aperture photometry 
was performed using the IDL APER procedure on each ob- 
ject. An aperture radius of L5" was adopted for each source 
(w 95% enclosed energy for soft band, w 90% enclosed en- 
ergy for hard band; aperture corrections were applied) except 
Jl 109+3736 and J1530+2310, for which the aperture radius 
was 3.0" because of their large numbers of detected X-ray 
counts (see Tabled no aperture corrections were applied to 
these two sources). The background region for each source 
was defined as an annulus with inner and outer radii of twice 
and three times the aperture radius. All background regions 
are free of X-ray sources. The upper limits upon X-ray counts 
for J0945-I-1009 were determined using the method of Kraft et 
al. (1991) at 95% confidence. Table|4]lists the X-ray counts 
in the three bands, as well as the band ratio (defined as the 
ratio between hard-band counts and soft-band counts) and ef- 
fective power-law photon index for each source. The effective 
power-law photon index was determined from the band ratio 
using the Chandra PIMM^3 tool, under the assumption of a 
power-law model with Galactic absorption only. 

Four archival sources (J1139-0201, J1604+4326, 
J2115+0001, and J2324+1443) were observed by Chandra 
in Cycles 5, 7, and 10. All of the quasars were detected 
by Chandra except J2115+0001. Similar Chandra data- 
reduction and processing procedures were performed for 
these objects. Three of them (J1139-0201, J2115+0001, 
and J2324+1443) were targeted in their Chandra obser- 
vations, for which the aperture radius was set to be 1.5". 
J1604+4326 was serendipitously covered by the ACIS-I 
detector in two Chandra observations. We measured the 
X-ray counts individually for these two observations, and 
then calculated the mean count rate and flux in the soft band. 
This source did not show significant variability (< 12%) 
between its two Chandra observations (w 10 hours apart 
in the quasar rest frame). The aperture radius (2.6") for 
this source was determined to be the 95% enclosed-energy 
radius at 1.497 keV based on the point spread function 
(PSF) of the ACIS detector at an off-axis angle of 3.7'. The 
Chandra observations of J21 15+0001 and J2324+1443 were 
briefly reported in Plotkin et al. (2010b), and our results are 
consistent with theirs. 

One archival source (J1013+4927) was serendipitously 
covered by XMM-Newton on 2004 April 23. Data reduction 
and processing were performed using standard XMM-Newton 
Science Analysis System (vlO.0.0) routines. This source 
is undetected by both the MOS and pn detectors using the 
EBOXDETECT procedure. Visual inspection of the images ver- 
ifies the non-detection of this source. We only used the data 
from the MOS detectors because they have higher angular reso- 
lution, which enables more reliable count extraction and back- 
ground estimation. The events files were filtered by removing 
background flaring periods (12% of the total exposure time) in 
which the count rate exceeded 0.35 s^' for events with ener- 
gies above 10 keV. The aperture for photometry (49.9" radius) 
was taken to be the 90% enclosed-energy radius at 1.5 keV 
based on the PSF of the MOS detectors at an off-axis angle of 
6.4'. The upper limits upon X-ray counts were determined to 
be 3\/N, where is the total counts within the aperture. 

Table|5]lists the key X-ray, optical, and radio properties of 

' ^ http ://cxc . harvard. edu/toolkit/ pimms .j sp 
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our low-redshift WLQs: 

Column (1): The SDSS equatorial coordinates (J2000) for 
the source. 

Column (2): The apparent /-band magnitude of the source 
using the SDSS quasar catalog BEST photometry. 

Column (3): The absolute /-band magnitude for the source. 
Mi, from the SDSS DR7 quasar catalog (Schneider et al. 
2010), calculated by correcting for Galactic extinction and 
assuming a power-law spectral index of ay = —0.5 (e.g., 
Vanden Berk et al. 2001). 

Column (4): The Galactic neutral hydrogen column density 
in units of 10^" cm"^, obtained with the Chandra COLDE?^3 
tool. 

Column (5): The count rate in the observed-frame soft 

http://cxc.harvard.edu/toolkit/colden.jsp 



X-ray band (0.5-2.0 keV), in units of 10 ^ s ^ For the two 
off-axis sources (J1013H-4927 and J1604H-4326), the count 
rate (or upper limit) is corrected for vignetting using exposure 
maps. 

Column (6): The Galactic absorption-corrected flux 
in the observed-frame soft X-ray band in units of 
10"'"^ erg cm"2 s"', obtained with the Chandra PIMMS tool. 
An absorbed power-law model was utilized with a photon in- 
dex r = 2, which is typical for quasars, and the Galactic neu- 
tral hydrogen column density for each source {Nh, given in 
Column 4). 

Column (7): The Galactic absorption-corrected flux density 
at rest-frame 2 keV in units of 10^^^ erg cm^^ s^' Hz^*, 
obtained with the Chandra PIMMS tool. 

Column (8): The logarithm of the quasar X-ray luminosity 
in the rest-frame 2-10 keV band corrected for Galactic ab- 
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sorption. 
Column (9): 



2500 A in units of 10^27 gj.g^.jj^-2g-i jj^-i 



The continuum flux density at rest-frame 

, from the SDSS 
quasar spectral property catalog in Shen et al. (2011). 

Column (10): The logarithm of the monochromatic lumi- 
nosity at rest-frame 2500 A, derived from the flux density at 
rest-frame 2500 A. A cosmological bandpass correction is uti- 
lized. 

Column (11): The X-ray-to-optical power-law slope, given 
by 

^ log(/2kev//2500 A) ^ q 334 j^g ( A}^] 



log(V2kev/V2500 a) 



/2500 A 



The flux density is measured per unit frequency. 

Column (12): Aap x, a parameter assessing the relative 
X-ray brightness (see defined as 



AOJqx ^ox(measured) ^ox(expected) • 



(2) 



The expected oCqx for a typical radio-quiet quasar is calculated 
using the 0Cox-i2500 A correlation given as Equation (3) of Just 
et al. (2007). The statistical significance of AOox (given in 
parentheses) is in units of a, which is obtained from Table 5 
of Steffen et al. (2006) as the RMS for Oqx of quasars with 
several ranges of luminosity. 

Column (13): The factor of X-ray weakness, derived from 
the Attox values in Column (12), quantifying the X-ray weak- 
ness of our sources compared to a typical radio-quiet quasar 
with similar UV/optical luminosity, calculated as /x-weak = 
jQ-Aa„x/o.384 ^ 403-Aaox A source with Aaox= -0.384 has 
an X-ray flux only sa 10% that of typical quasars, correspond- 
ing to an X-ray weakness factor of w 10. 

Column (14): The optical-to-radio power-law slope, given 
by 

log(/5GHz//2500 a) ^t,) 



Oro ^ 



log(V5GHz/V25ooA)' 



The values of /2500 A given in Column (9). The values of 
/s GHz were calculated using a radio power-law slope of Uy = 
-0.8 (e.g., Falcke et al. 1996; Barvainis et al. 2005) and a 
flux at an observed-frame wavelength of 20 cm, /20 cm- For 
sources detected by the FIRST survey (Faint Images of the 
Radio Sky at Twenty-Centimeters; Becker et al. 1995), /20 cm 
was taken from the FIRST source catalog. For sources cov- 
ered but not detected by the FIRST survey, the upper limits for 
radio flux density were placed as /20 cm < 0.25 + (5c7nns) mJy, 
where Cnns is the RMS noise of the FIRST survey at the ob- 
ject's coordinates (see §5.3.1 of Plotkin et al. 2010a for more 
details). For the two sources not covered by the FIRST sur- 
vey (J21 15H-0001 and J2324H- 1443), the radio flux density was 
measured via targeted VLA observations (see §3.1 of Plotkin 
et al. 2010b). 

The Oro parameter is related to the commonly used radio- 
loudness parameter, = /s ghz//44oo A (^-g-, Kellermann 
et al. 1989), by the following equation: 



«10 



log[iJ(2500/4400) 

log(V5GHz/V2: 



-0.186 log/? -0.023, (4) 



500 A) 



where we use 



-0.5. Therefore, we have aio> —0.21 



for radio-quiet quasars (R < 10); -0.39 < aro< -0.21 for 
radio-intermediate quasars (10 < R < 100); and 0!ro< —0.39 
for radio-loud quasars (R > 100). All the low-redshift WLQ 
candidates in our sample are radio quiet. Note that the radio 



and optical observations of our sources are non-simultaneous. 
While the non-simultaneity does not generally alter the clas- 
sification of typical radio-quiet quasars, it may significantly 
affect that of BL Lac objects because of their rapid and large- 
amplitude variability. 

5. RESULTS AND DISCUSSION 
5.1. Relative X- ray B righ tness 

The Attox parameter (see Eqn. |2] for definition) is utilized 
to assess the X-ray brightness of a quasar relative to typical 
radio-quiet quasars with similar UV luminosity. We compare 
the Attox distribution of our low-redshift, radio-quiet WLQ 
candidates (see Fig.|6]l to that of the 132 radio-quiet, non-BAL 
quasars in Sample B of Gibson et al. (2008)0 which repre- 
sent typical radio-quiet SDSS quasars. All of the 132 Sam- 
ple B quasars are X-ray detected. The Peto-Prentice test (e.g., 
Latta 1981), implemented in the Astronomy Survival Analy- 
sis package (ASURV; e.g., Lavalley et al. 1992), is used to 
assess whether our low-redshift WLQ candidates follow the 
same Attox distribution as that for typical quasars (see results 
in Table |6l). We prefer the Peto-Prentice test to other pos- 
sible similar tests because it is the least affected by the fac- 
tors of different censoring patterns or unequal sizes of the two 
samples which exist in our case. We also compare the Attox 
distribution of high-redshift, radio-quiet WLQs in Shemmer 
et al. (2006, 2009) to that of our low-redshift, radio-quiet 
WLQ candidates and that of typical SDSS quasars (also see 
Fig. i] and Table©. 

The Attox distribution of our low-redshift, radio-quiet WLQ 
candidates is significantly different from that of typical SDSS 
quasars. The probability of null-hypothesis (two samples fol- 
lowing the same distribution) is only 6.3 x 10^^. This result is 
mainly due to the presence of a skew tail of X-ray weak WLQs 
(see Fig. |6]l. Seven out of the 1 1 objects in our sample of 
low-redshift, radio-quiet WLQ candidates have Attox < —0.2, 
giving a fraction of X- ray weak objects of (64+^^)% (68% 
confidence level). The mean Attox value for the low-redshift, 
radio-quiet WLQ candidates is —0.214 ± 0.078, calculated 
with the Kaplan-Meier estimator also implemented in the 
ASURV package, while that for the Sample B quasars is 
-0.001 ± 0.011. The Aaox distribution of the nine high- 
redshift, radio-quiet WLQs in Shemmer et al. (2006, 2009) is 
also different from that of typical radio-quiet SDSS quasars, 
but less significantly (the probability of null-hypothesis is 
7.1 X 10-^)0 Five of the nine objects in the high-redshift, 
radio-quiet WLQ sample have Attox < —0.2, giving a frac- 
tion of X-ray weak objects of (56^^4)% (68% confidence 
level); we note that the fraction could be somewhat higher 
(6/9) owing to the weak X-ray upper limit for J1237H-6301. 
The mean Attox value for the high-redshift, radio-quiet WLQs 
is —0.144 ± 0.075. As expected the combined low-redshift 

" We used an improved version of tlie Sample B quasars in Gibson et al. 

(2008) from which we further removed seven BAL quasars (see Footnote 16 
inWuetal. 2011). 

'** This result is somewhat inconsistent with the finding by Shemmer et al. 

(2009) that their high-redshift WLQs have a similar AOox distribution to that 
of typical SDSS quasars. Shemmer et al. (2009) used the Kolmogorov- 
Smirnov test and ignored the two high-redshift, radio-quiet WLQs with AOox 
upper limits. We include all nine high-redshift, radio-quiet WLQs since the 
Peto-Prentice test can properly treat censored data. The utilization of an im- 
proved Sample B (see Footnote 17) does not substantially contribute to the 
inconsistency here. The Peto-Prentice test using all nine high-redshift, radio- 
quiet WLQs and the original Sample B (as used by Shemmer et al. 2009) 
provides a null-hypothesis probability of 1.53 x 10^^. 
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Fig. 6. — Distribution of AOox values for the WLQ candidates, compared to that of the 132 radio-quiet, non-BAL quasars in Sample B of Gibson et al. (2008). 
The red histograms and red leftward arrows represent the low-redshift sources in our sample which are detected (8 sources) and undetected (3 sources) in X-rays, 
respectively. The blue histogram and blue leftward aiTows represent the high-redshift sources from Shemmer et al. (2006, 2009) which ai'e detected (7 sources) 
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1811 is shown as a red asterisk. The dashed vertical line shows Ao!ox= 0. Note the many (> 50%) WLQs with AOa^ < —0.2. 



and high-redshift, radio-quiet WLQ sample (mean Ao£ox value 
of —0.187 ± 0.056) also follows a different AoCqx distribution 
from that of typical SDSS quasars (null-hypothesis proba- 
bihty of 4.5 x 10"^). The Attox distribution of low-redshift, 
radio-quiet WLQ candidates is consistent with that of high- 
redshift, radio-quiet WLQs (null-hypothesis probability of 
0.35), though the sample sizes being compared are limited. 

5.2. Classifying Radio-Quiet WLQs 

To investigate the multi-band properties of low-redshift, 
radio-quiet WLQ candidates, we plotted the sources of our 
sample in an Oro-Olox diagram (Fig. |7]) along with the high- 
redshift, radio-quiet WLQs in Shemmer et al. (2006, 2009), 
the BL Lac sample in Shemmer et al. (2009), and the Sam- 
ple B quasar in Gibson et al. (2008). The low-redshift, radio- 
quiet WLQ candidates have similar multi-band properties to 
those of high-redshift, radio-quiet WLQs. They are generally 
much fainter in radio and X-rays than most of the BL Lac 
objects. The weak emission lines of low-redshift, radio-quiet 
WLQ candidates are therefore not likely due to the dilution 
by relativistically boosted continua as for BL Lac objects (see 
discussion in §4.1 of Shemmer et al. 2009). However, it is 
possible that a small percentage of the WLQ candidates ac- 
tually belong to the radio-faint tail of the BL Lac population 
(see below). 

The population of radio-quiet WLQ candidates (both at low 



redshift and high redshift) has a wide dispersion of relative 
X-ray brightness and UV emission-line properties. Motivated 
by their Attox distribution, their emission-line properties dis- 
cussed below, and observations of related objects (e.g., Wu 
et al. 201 1), we will discuss them in three groups. 

The majority of WLQ candidates are not only X-ray weaker 
than BL Lac objects, but also weaker than typical radio- 
quiet SDSS quasars. These WLQ candidates may belong to 
the notable class of X-ray weak quasars termed "PHL 1811 
analogs" which were recently studied in detail by Wu et al. 
(2011). The PHL 1811 analogs generally have weak and 
highly blueshifted high-ionization lines (e.g., C IV, Si iv), 
weak semi-forbidden lines (e.g., C III]), and strong UV Fe II 
and/or Fe III emission. Some of our low-redshift, radio-quiet 
WLQs have similar UV emission-line properties to those of 
PHL 1811, as listed below: 

1. J0812H-5225 (Aaox= -0.42) has weak C III] and strong 
Fe II emission. 

2. J0945H-1009 (Aaox < -0.34) has weak C iv and C ill] 
emission lines. Its C IV line is highly blueshifted (« 
-7000kms-i). 

3. J1252H-2640 (Aaox= -0.39) has weak C III] and strong 
Fe II emission. 

4. Jl 139-0201 (Aaox= -0.38) has weak and highly 
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blueshifted (sa -2950 km s"') C IV emission, weak 
C III] emission, and strong Fe III emission. 

A high-redshift, radio-quiet WLQ J1302H-0030 (Aaox= 
—0.38) also has a weak and highly blueshifted C IV emis- 
sion line (DS09; Wu et al. 2011). All of the above men- 
tioned sources are X-ray weak by a factor of > 7 (see Ta- 
ble [Sj. Fig. |8] shows the distribution of our WLQ candi- 
dates in the Attox-C IV blueshift- REW(C IV) parameter 
space. J0945H-1009, Jl 139-0201 and J1302H-0030 are sim- 
ilar to PHL 1811 analogs in this diagram. Based on the model 
in §4.6 of Wu et al. (2011), these PHL 1811 analogs may 
have high-ionization shielding gas with large column density 
and a large covering factor of the BELR which blocks most of 
the ionizing photons, resulting in weak high-ionization emis- 
sion lines. If a quasar of this kind is viewed through the 
BELR and shielding gas, it would be an X-ray weak WLQ 
with weak and highly blueshifted high-ionization lines (e.g., 
C IV). Based on the estimate in §4.6 of Wu et al. (2011), 
PHL 1811 analogs should make up « 30% of the total WLQ 
population. However, our sample of low-redshift, radio-quiet 
WLQ candidates appears to have a higher fraction (> 50%) of 
PHL 1811 analogs, which may indicate our sample has some 
selection bias toward X-ray weak WLQ candidates. This bias 
could perhaps be the result of a more strict criterion upon 
the strengths of emission lines for most sources in our sam- 
ple (REW < 5 A). Quasars with weaker emission lines (e.g.. 



C IV) are perhaps more likely to be weak in X-rays (e.g., see 
§4.5 of Wu et al. 2011). The apparently higher fraction of 
PHL 1811 analogs in our sample than that in Wu et al. (201 1) 
may perhaps also be caused simply by small-sample statistics. 
A Fisher's exact test (Fisher 1922) gives an 11.1% probability 
for the different fractions of PHL 1811 analogs among these 
two samples under the null hypothesis (i.e., the two samples 
have the same fraction of PHL 1811 analogs). 

Some of our WLQ candidates have similar X-ray brightness 
to that of typical radio-quiet quasars (—0.2 < AOox ^ 0.2). 
Their high-ionization lines are also weak, but perhaps not 
highly blueshifted (e.g., see Fig. [8] for J1408H-0205). Some 
of them (e.g., J1612-I-51 18) have very weak UV Fe II and/or 
Fe III emission. J1612-I-5118 does seem to have a highly 
blueshifted C IV line, for which the reason is unclear. How- 
ever, the C IV line of this source is close to the blue border 
of its SDSS spectral coverage. Further UV spectroscopy with 
better C IV coverage is needed to confirm its C IV blueshift. 
Based on the model in Wu et al. (201 1), these sources are sim- 
ilar to PHL 1811 analogs physically, but they are viewed at 
different orientations. These sources are observed along lines 
of sight that avoid the shielding gas and the BELR. Therefore 
they appear normal in X-rays. Their high-ionization lines are 
generally not highly blueshifted. 

In our WLQ candidate sample, one source (Jl 109-1-3736) is 
remarkably strong in X-rays. It also shows similar UV/optical 
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Fig. 8. — REW(C IV) plotted against the C IV blueshift for our radio-quiet WLQ candidates (filled squares for X-ray detected sources, filled upside-down 
triangles for X-ray undetected sources), PHL 1811 (asterisk), radio-quiet PHL 1811 analogs in Wu et al. (2011) (stars), and radio-quiet, non-BAL quasars in 
Sample B of Gibson et al. (2008a) (circles). These sources are c olor-coded according to their Attox values (three color bins are used, corresponding to the X-ray 
weak. X-ray normal and X-ray strong sources described in ^5.2| respectively). The color bar shows the Aa^x range for each color. Source names for WLQs are 
labeled in the format of 'Ihhiwif for brevity. The grey dots show the 13,582 radio-quiet quasars in Sample A of Richards et al. (201 1; see their Fig. 7). 



spectral properties to those of BL Lac objects. This source 
may belong to the radio-faint tail of t he BL Lac popula- 
tion; we will discuss it further in ^5.5\ It is worth noting 
that the division of our radio-quiet WLQ candidates into the 
three groups discussed above (as shown in Fig. |7]i is some- 
what arbitrary. We do have some "border-line" sources with 
AOox ~ ±0.2 (e.g., J1212H-5341). It is difficuh to classify 
these sources clearly based on current information. 

5.3. The Infrared-to-X-ray SEDs of the Radio-Quiet WLQ 
Candidates 

For the purpose of investigating further the multi- 
wavelength SEDs of our low-redshift, radio-quiet WLQs, 
we gathered photometry for our sample from the following 
bands; (1) near- and mid-infrared from WISE (The Wide- 
field Infrared Survey Explorer; Wright et al. 2010); (2) near- 
infrared from 2MASS (The Two Micron All Sky Survey; 
Skrutskie et al. 2006); (3) optical from the SDSS; (4) UV 
from GALEX (The Galaxy Evolution Explorer; Martin et al. 
2005); and (5) X-rays from this work. Fig.|9]shows the SEDs 



of the five low-redshift WLQs in our sample which have the 
best multi-band coverage. We also include the SED of the 
radio-quiet BL Lac candidate J1109-I-3736 (see more discus- 
sion in ^3^. A key point to keep in mind is that these multi- 
band observations are non-simultaneous. The SEDs are there- 
fore subject to potential distortions due to variability. 

We examined the WISE and GALEX image tiles by eye 
to identify potential cases of source blending, confusion, or 
incorrect matching caused by the low angular resolution of 
WISE and GALEX. None of the sources in our sample is sub- 
ject to these kinds of problems. The 2MASS magnitudes 
in the SDSS DR7 quasar catalog were utilized; this cata- 
log provided aperture photometry for additional sources de- 
tected down to 2(7 (see §5 of Schneider et al. 2010). For 
the sources without detections at > 2(7, we adopted flux up- 
per limits obtained following the same photometry procedure 
(C. M. Krawczyk & G. T. Richards 2011, private communi- 
cation). The first five sources in Fig. |9] are from the groups 
of X-ray weak and X-ray normal WLQ candidates discussed 
above. Their mid-infrared-to-UV SEDs are generally con- 
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sistent with the composite SEDs of typical SDSS quasars in 
Richards et al. (2006), and they are significantly different 
from the SEDs of BL Lac objects (see the dotted parabolic 
line in the top-left panel of Fig.|9l e.g., Nieppola et al. 2006). 
We also investigated the SEDs for the WISE-coveied radio- 
quiet objects cataloged in Plotkin et al. (2010a) which do 
not have sensitive X-ray coverage (see the Appendix). The 
majority of them also have mid-infrared-to-UV SEDs con- 
sistent with those of typical radio-quiet quasars. Lane et al. 
(201 1) obtained similar results for their high-redshift WLQs; 
the composite SED of their high-redshift WLQs is inconsis- 
tent with SEDs of BL Lac objects. For one source in our 
sample, J0945H-1009, the flux in the UV band is lower than 
for typical SDSS quasars (see the GALEX data points in the 
top-right panel of Fig. |9]l. The UV deficiency of this source 
may be caused by Lya-forest intervening absorption. How- 
ever, Laor & Davis (2011) argued that such intervening ab- 
sorption is not significant (~ 1 1 % at most) for this source with 
z = 1.66. The near-infrared-to-UV SED of J0945H-1009 can 
be well fitted with their local black-body model for a cold ac- 
cretion disk. 

5.4. X-ray Spectral Properties of Low-Redshift WLQ 
Candidates 

Most of the low-redshift, radio-quiet WLQ candidates do 
not have sufficient X-ray counts for an individual X-ray spec- 
tral analysis. We therefore investigate the average X-ray spec- 
tral properties for these sources via stacking analyses and joint 
fitting. 

A stacked spectral analysis was performed for the six 
low-redshift, radio-quiet WLQ candidates with Attox < 
-0.3 (J0812H-5225, J0945H-1009, Jl 139-0201, J1252H-2640, 
J2115H-0001, and J2324H-1443). These sources are the weak- 
est in X-rays among the full sample of low-redshift WLQ can- 
didates. The detected sources have similar numbers of X-ray 
counts, so that any one of them will not dominate the stacking 
analysis. These six sources span a relatively wide range of 
redshift (z = 1.15-2.50), and thus the observed-frame bands 
of each source correspond to different energy ranges in the 
rest frame. However, under the assumption of a simple power- 
law spectral model, one can stack the X-ray counts to ob- 
tain the average effective power-law photon index. We added 
the X-ray counts of these sources in the observed-frame soft 
band and hard band, respectively. The numbers of total net 
counts are 14.4^3 ^ in the soft band and 5.2^2 2 the hard 
band (68% confidence level), and the resulting band ratio is 
0.36to fs- With the average Galactic neutral hydrogen col- 
umn density of these sources (A^h = 3.50 x 10^" cm^^), the 
band ratio was converted to an effective power-law photon in- 
dex r = 1.66^05 J. The average X-ray spectrum of the X-ray 
weak low-redshift WLQs is perhaps somewhat harder than 
that for typical radio-quiet quasars (F w 2), but it is consistent 
within the error bars. This average X-ray spectrum is likely 
softer than that of the PHL 1811 analogs at z = 2.19-2.38 
(F = 1.10;S3;^f,) in Wu et al. (2011) which was also obtained 
via a stacking analysis, but consistent within 2a. Both stack- 
ing analyses suffer from large uncertainty due to limited X-ray 
counts. For a sample combining the six X-ray weak WLQs 
analyzed here (which are likely to be PHL 1811 analogs) and 
the radio-quiet PHL 1811 analogs in Wu et al. (2011), the 
average X-ray spectrum has a flat effective power-law photon 
index of F — 1 .35^^ 3;. Deeper X-ray observations are neces- 
sary to give tighter constraints on the X-ray spectral properties 



of these X-ray weak WLQ candidates. 

Two sources (J0945H-1009 and J21 15H-0001) are undetected 
by Chandra. Adding the X-ray counts of these two sources 
cannot generate a stacked source that would be detected by 
Chandra, because both of these two sources have zero X-ray 
counts. However, we are able to obtain a tighter average con- 
straint on their X-ray brightness via stacking analysis. The 
upper limit upon the soft-band count rate of the stacked source 
is 2.51 X lO^"* s^^ Average values of redshift. Galactic A^h- 
and /2500 A ^re adopted in the following calculation. The up- 
per limit upon the average flux density at rest-frame 2 keV 
is 5.16 X lO^-'^ erg cm^^ s^' Hz^' under the assumption of 
the Galactic-absorbed power law with F = 2. The upper lim- 
its upon aox and Attox are calculated to be Oox < —2.21 and 
AOox < —0.50. Therefore, these two sources are X-ray weak 
by a factor of > 20 on average. 

For the two X-ray normal, low-redshift WLQs 
(J1604H-4326 and J1612H-51 18), we performed joint fit- 
ting to study the average X-ray spectral properties of these 
sources. The X-ray spectra were extracted from apertures of 
3" radius centered on the X-ray positions of these sources 
via the standard ClAO routine PSEXTRACT. The background 
spectra were extracted from annular regions with inner radii 
of 6" and outer radii of 9", which are free of X-ray sources. 
Two spectra were extracted individually for J1604H-4326 from 
its two Chandra observations. Spectral fitting was performed 
with XSPEC V12.6.0 (Arnaud 1996). The C-statistic (Cash 
1979) was used in the spectral fitting instead of the standard 
statistic because the C-statistic is well suited to the limited 
X-ray counts in our analysis (e.g., Nousek & Shue 1989). 
We fit the spectra jointly using a power-law model with a 
Galactic absorption component represented by the wabs 
model (Morrison & McCammon 1983). We also used another 
model similar to the first, but adding an intrinsic (redshifted) 
neutral absorption component, represented by the zwabs 
model. Both sources were assigned their own values of 
redshift and Galactic neutral hydrogen column density; the 
Galactic column density was fixed to the values calculated 
with GOLDEN (Column 4 of Table |5]l. The joint fitting 
results are shown in Table [7] The quoted errors or upper 
limits are at the 90% confidence level for one parameter of 
interest (AC = 2.71; Avni 1976; Cash 1979). The average 
X-ray spectral properties of these two sources are similar to 
those of typical radio-quiet quasars. They have an average 
photon index (F = 2.07^q 3q), consistent with that of typical 
radio-quiet quasars (F sa 2). The average photon index is also 
consistent with those from their band-ratio analyses. We did 
not find evidence of strong intrinsic neutral absorption for 
J1604H-4326 and J1612H-5118 (A^h < 1-58 x 10^2 cm^^); 
the spectral fitting quality was not improved after adding the 
intrinsic neutral absorption component. 

Two sources (Jl 109H-3736 and J1530H-2310) have sufficient 
X-ray counts for individual spectral analysis. We will dis- 
cuss the X-ray spectral properties of J1530H-2310 here and 
leave the radio-quiet BL Lac candidate J1109H-3736 for the 
next subsection. The X-ray spectrum of J1530H-2310 was ex- 
tracted following the same procedure as for J1604H-4326 and 
J1612H-5118 (see above). The spectrum of J1530H-2310 was 
first grouped to have at least 10 counts per bin (see Fig. [TOb . 
We used the same spectral models as those for the joint fit- 
ting above. The standard statistic was utilized. The fit- 
ting results are also shown in Table |7] Fig. [TO] shows the 
X-ray spectrum, the best-fit power-law model with Galactic 
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Fig. 10.— X-ray spectra of SDSS J 1530+23 10 (top panel; binned to a min- 
imum of 10 counts per bin) and SDSS Jl 109+3736 (bottom panel; binned to 
a minimum of 20 counts per bin) fitted with a power-law model with Galactic 
absorption. The residuals are shown in units of a. The inset of each panel 
shows contours of the photon index vs. intrinsic neutral hydrogen column 
density parameter space, at confidence levels of 68%, 90%, and 99%, respec- 
tively. 

absorption, and the contour plot of the F — A^h parameter space 
for the spectral model with intrinsic neutral absorption for 
J1530H-2310. The X-ray spectral properties of J1530H-2310 
are similar to those of J1604H-4326 and J1612H-5118. The 
photon index of J1530H-2310 from the spectral fitting (F = 
2.11^034) is consistent with that from the band-ratio anal- 
ysis (see Table nil. J1530H-2310 also shows no evidence of 
strong intrinsic neutral absorption (A^h ^ 2.67 x 10^^ cm^^). 
In summary, the X-ray spectral properties of X-ray normal 
low-redshift WLQs are consistent with those of typical radio- 
quiet quasars. Shemmer et al. (2009) found similar results 
for their high-redshift, radio-quiet WLQs. Their sources have 
a somewhat harder average X-ray spectrum (F = IMtolj), 
perhaps because they fit both X-ray weak and X-ray normal 
WLQs jointly. Shemmer et al. (2009, 2010) also performed 
individual X-ray spectral analysis on two high-redshift, radio- 
intermediate WLQs (J1141H-0219 and J1231H-0138). 

5.5. Jl 109+3736: The Radio-quiet BL Lac Candidate 

Jl 109H-3736 is a radio-quiet BL Lac candidate based on its 
multi-band properties. It is strong in X-rays by a factor of 
w 6.3 (Oox = —1-10, Attox = 0.31). Its oCox value is similar 
to that of the majority of the BL Lac population (see Fig.|7ll. 
X-ray spectral analysis was carried out for J1109H-3736 fol- 
lowing similar procedures as those in i)5.4| The Chandra 
spectrum of J1109H-3736 was grouped to have at least 20 
counts per bin (see Fig. [TOl i. The best-fit parameters are 
listed in Table |7] The X-ray spectrum, the best-fit power- 
law model with Galactic absorption, and F — A^h contours for 
J1109H-3736 are also shown in Fig. [TO] The best-fit photon 



index of Jl 109H-3736 (F = 1 .77 ±0.14) is consistent with that 
from the band-ratio analysis. This photon index value indi- 
cates a harder X-ray spectrum than those for the majority of 
the high-energy peaked BL Lac objects (HBL){3 which have 
a mean photon index F « 2.2, but it is still consistent with 
the broad distribution of HBL photon-index values (e.g., see 
the bottom panel of Fig. 1 in Donato et al. 2005). Although 
this source has a high X-ray count rate, we do not expect 
strong photon pile-up effects because a 1/2 subarray mode 
was used for its Chandra observation. There is no evidence 
for significant intrinsic absorption in the J1109H-3736 spec- 
trum (A^H ^ 1-4 X 10^^ cm"^). 

The SDSS mectrum of Jl 109H-3736 (see the bottom-right 
panel of Fig. |2]i shows a strong power-law continuum without 
any detectable emission lines, in particular no Balmer lines. 
Note that some of the X-ray weak, low-redshift WLQ candi- 
dates (e.g., PHL 1811; see Leighly et al. 2007b) have fairly 
strong Balmer lines. The strength of its Ca II H/K breatPI 
(C = 0.185) indicates a relatively small contribution to the 
SDSS spectrum from the host galaxy. J1109H-3736 shows 
moderate X-ray variability between its Chandra observation 
and the epoch of the RO SAT AW Sky Survey in 1990 Novem- 
ber (RASS; Voges et al. 1999). It was not detected by RASS. 
The upper limit for Oox is estimated to be ttox < — 1 -29, show- 
ing a factor of > 3 variation compared to its Chandra obser- 
vation. The quasar should have been detected in the RASS if 
it had the same X-ray brightness as that in the Chandra epoch 
(expecting w 20 counts in a 250 s RASS observation). There 
is also evidence for optical variability of J1109H-3736. This 
source was optically identified during the Cambridge APM 
(Automated Plate Measuring Machine) scans of POSS (Palo- 
mar Observatory Sky Survey) plates (e.g., McMahon et al. 
2002). We converted the APM magnitudes in the O and 
E bands {O ~ 20.76, E = 19.12) into a B-band magnitude 
(B = 20.57) using Equations (1-3) in McMahon et al. (2002). 
We also converted the SDSS magnitudes to a Z?-band magni- 
tude (B ~ 18.88) using the equations in Table 1 of Jester et al. 
(2005). The B-band magnitude difference is 1.69, indicating a 
factor of 5 variability over a time span of ^ 50 years. This 
source shows greater variability than most quasars during this 
time span (e.g., see Fig. 24 of Sesar et al. 2006). Therefore, 
J1109H-3736 is a BL Lac candidate on the radio-faint tail of 
the full BL Lac population. 

The broad-band SED of Jl 109H-3736 is shown in Fig. |9] It 
is relatively weak in the near infrared compared to its SDSS 
brightness. However, the SED profile of this source could 
be strongly affected by its variability (see the comparison be- 
tween its SDSS and POSS fluxes in Fig.|9]) given its possible 
BL Lac nature. Infrared photometry from WISE, which will 
be available in the WISE full data release^ will be helpful 
because the infrared-to-UV SEDs of BL Lac objects are very 
different from those of typical quasars (e.g.. Lane et al. 201 1). 

" J 11 09+3736 would be classified as an HBL if it were indeed a 
BL Lac object, based on the criteiion < 0.75 of Padovani & 

Giommi (1995), where a^^ is the radio-to-X-ray spectral index, 0!fx = 
-0.130 log(/i kev//5 GHz). The Cfrx value of Jl 109+3736 is estimated to 
be 0.46. 

The strength of the Ca 11 H/K break is defined as the fractional change 
of the average flux densities in continuum regions blueward (3750-3950 A) 
and redward (4050-4250 A) of the H/K break (see Equation 1 of Plotkin et al. 
2010a). 

^' Jl 109+3736 is not covered by the currently avail- 

able WISE preliminary data release. See coverage map at 

|http://wise2.ipac. caltech. edu/docs/release/prelim/| 
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Future polarization measurements of Jl 109+3736 will also 
be useful in understanding its nature, since BL Lac objects 
are usually highly polarized in the optical/UV band. How- 
ever, polarimetry surveys of optically selected BL Lac sam- 
ples (e.g.. Smith et al. 2007; Heidt & Nilsson 2011) did not 
find any highly polarized radio-quiet BL Lac candidates, in- 
dicating this kind of source, if it indeed exists, should make 
up only a small fraction of the total population of radio-quiet 
WLQ candidates. 

6. SUMMARY AND FUTURE STUDIES 

We have compiled a sample of 1 1 radio-quiet WLQ can- 
didates with z = 0.4-2.5 and presented their X-ray and mul- 
tiwavelength properties. These sources are mainly selected 
from the catalog of radio-quiet, weak-featured SDSS quasars 
in Plotkin et al. (2010a). Six of them were observed in 
new Chandra Cycle 12 observations, while five have archival 
Chandra or XMM-Newton coverage. Our main results are 
summarized as follows: 

1. All newly observed low-redshift, radio-quiet 
WLQ candidates are detected by Chandra, ex- 
cept for J0945+1009. Three sources (JOS 12+5225, 
J0945+1009, and J1252+2640) are X-ray weak by fac- 
tors of > 8-12 compared to typical quasars with similar 
optical/UV luminosity. Two sources (J1530+2310 and 
J1612+5118) are X-ray normal, while the other one 
(Jl 109+3736) is X-ray strong by a factor of 6.3. See 
SI 

2. Three (Jl 139-0201, J1604+4326, and J2324+1443) 
of the five low-redshift, radio-quiet WLQ candidates 
with sensitive archival X-ray coverage are detected in 
X-rays, while two (J 101 3+4927 and J2 115+0001) do 
not have X-ray detections. All of the five archival 
sources are X-ray weak by factors of > 3-12. See SI 

3. The distribution of relative X-ray brightness for our 
low-redshift, radio-quiet WLQ candidates is signifi- 
cantly different from that of typical radio-quiet quasars. 
Our sample has a highly statistically significant excess 
of X-ray weak sources. About 64% (7/1 1) of the low- 
redshift, radio-quiet WLQs and about 56% (5/9) of the 
high-redshift, radio-quiet WLQs are X-ray weak. The 
X-ray weakness that is commonly found within WLQ 
samples may well be the driver of the weak broad-line 
emission. Therefore, X-ray weakness provides an im- 
portant clue for understanding the nature of WLQs. See 

a 

4. The X-ray weak sources (Aaox< —0.2) in our low- 
redshift WLQ sample are likely to be PHL 1811 
analogs (see Wu et al. 201 1). Some of them show sim- 
ilar UV emission-line properties to those of PHL 1811 
(weak and highly blueshifted high-ionization lines, 
weak semi-forbidden lines, and strong UV Fe emis- 
sion). A stacking analysis of these sources indicates 
the average effective power-law photon index to be 
r = IMtlfi (68% confidence level). See ^ and 

5. Sources with —0.2 < Attox 0.2 have multi-band prop- 
erties that suggest they are X-ray normal WLQs. They 
also have weak high-ionization lines, while some of 
them do not have highly blueshifted high-ionization 



lines. Some X-ray normal sources in our sample have 
weak UV Fe emission. The average X-ray spectral 
properties of these sources are similar to those of typi- 
cal SDSS quasars. According to the unification model 
in §4.6 of Wu et al. (2011), these sources may have a 
similar geometry and physical nature to the P HL 1 8 1 1 
analogs, but with different viewing angles. See ii5.2l and 

6. The mid-infrared-to-UV SEDs of the X-ray weak and 
X-ray normal low-redshift WLQ candidates are gen- 
erally consistent with the composite SEDs of typical 
SDSS quasars, suggesting that these sources are not 
likely to be BL Lac objects with relativist ically boosted 
continua and diluted emission lines. See ^15.31 

7. The X-ray strong source Jl 109+3736 (Aaox> 0.3) may 
be a radio-quiet BL Lac object. It has similar X-ray 
brightness (aox= — 1 . 10) to those of typical BL Lac ob- 
jects. The SDSS spectrum of Jl 109+3736 shows a sim- 
ilar continuum to those of BL Lac objects. It also has 
shown moderate X-ray (a factor of > 3) and optical (a 
factor of « 5) variability in comparisons with archival 
data. A Chandra spectral analysis using a Galactic- 
absorbed power-law model gives F = 1 .77 ± 0. 14; thus 
its X-ray spectrum is harder than those of the majority 
of high-energy peaked BL Lac objects. There is no ev- 
idence of significant intrinsic X-ray absorption for this 
source. See ^5.5\ 

Future studies of larger samples of radio-quiet WLQ candi- 
dates will be helpful to clarify their nature. UV spectroscopy 
covering the Lya + N V and C IV regions of low-redshift, 
radio-quiet objects is necessary to study the REW distribu- 
tions of these two lines, and their REW correlations. This 
will provide insights toward a possible universal definition 
for WLQs at all redshifts, which should enable systematic 
studies of larger, unbiased, and more complete samples of 
WLQ candidates. For current WLQ studies, it is difficult 
to perform reliable mid-infrared-to-X-ray SED and correla- 
tion analyses because of the limited sample size. Further ac- 
cumulation of high-quality X-ray data will substantially en- 
large the sample sizes for such analyses. Deeper X-ray ob- 
servations are required to convert the X-ray flux upper lim- 
its into detections and thus to study the true overall distribu- 
tion of relative X-ray brightness for radio-quiet WLQ can- 
didates. High-quality X-ray spectroscopy should be able to 
reveal any X-ray absorbers in X-ray weak WLQ candidates, 
clarifying the cause of their X-ray weakness and the geometry 
of these quasars. More accurate measurements of the photon 
indices of their hard X-ray power-law spectra can better con- 
strain the L/L-Edd values of WLQ candidates (e.g., Shemmer 
et al. 2008) to see whether extremely high/low L/L-eaa is the 
physical cause of their weak broad emission lines. Therefore, 
radio-quiet WLQ candidates will be excellent targets of future 
missions with much higher X-ray spectroscopic capability, 
e.g., the Advanced Telescope for High Energy Astrophysics 
(A77/£7VA)0Near-infrared spectroscopy covering the H)3 re- 
gion of the low-redshift WLQ candidates will also allow esti- 
mation of their L/Lgdci values (e.g., Shemmer et al. 2010). 

Although most WLQ candidates likely do not have rela- 
tivistically boosted continua, this work suggests that a minor- 
ity of them may belong to the radio-faint tail of the BL Lac 

http://www.mpe.mpg.de/athena/home.php 
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population. Further growth of the high-quality multiwave- 
length database, especially in the infrared band, is crucial to 
study the broad-band SEDs of WLQ candidate s, wh ich could 
distinguish BL Lac objects from WLQs (see ^3^. Full re- 
lease of the WISE source catalog will greatly benefit SED 
studies of radio-quiet WLQ candidates. BL Lac objects often 
have large-amplitude variability and high polarization in the 
UV/optical band. Long-term UV/optical monitoring and po- 
larimetry of X-ray strong WLQ candidates will help to iden- 
tify their nature conclusively. 
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APPENDIX 

SPECTRAL ENERGY DISTRIBUTIONS OF ADDITIONAL RADIO-QUIET WLQ CANDIDATES 

The broad-band SEDs of WLQ candidates are able to provide useful insights into their nature. Lane et al. (2011) showed that 
the mid-infrared-to-UV SEDs of their high-redshift WLQs were consistent with those of typical quasars, but were significantly 
different from those of typical BL Lac objects. In f|53]we have obtained similar results for the X-ray weak and X-ray normal 
WLQ candidates in our low-redshift sample (see Fig.|9]l. However, it is possible to gain insights from SED measurements even 
when sensitive X-ray data are not available. Therefore, in this appendix we will further study the mid-infrared-to-UV SEDs of 
the radio-quiet weak-featured AGNs cataloged in Plotkin et al. (2010a) that do not have sensitive X-ray coverage. 

We cross-correlated the objects in Table 6 of Plotkin et al. (2010a) to the WISE source catalog in its preliminary data release, 
obtaining 30 VW^ii-detected sources (not including the four sources already discussed in ^3^. The WISE coverage of the 
remaining sources was checked with the online image tile look-up tool|3 An additional six sources with WISE coverage were 
identified. Among these six sources, two (J0901H-3846 and J1409— 0000) were identified as stars based on their SDSS spectra 
and proper-motion data; these two sources were removed from our SED study. For the remaining four sources, we examined their 
WISE image tiles. Two of them (J0755H-3525 and J1448H-2407) have WISE detections in all four bands that lie below the catalog 
limit. We performed aperture photometry (using a standard 8.25" aperture radius) and obtained their fluxes by scaling their 
counts in the aperture to those of nearby sources (within 60" separation) appearing in the WISE catalog. The other two sources 
(J0857H-2342 and J1541H-2631) were not detected in their WISE images. Following the standard WISE photometry procedure, 
we calculated their flux upper limits at a 95% confidence level by adding the aperture flux measurement plus two times the 
uncertainty. For J0857H-2342, we could only obtain a flux upper limit in the w3 band, because there are nearby bright sources in 
its aperture in the other band im ages. The photometric data in the other near-infrared-to-UV wavebands were obtained following 
the same methods described in g53\ The X-ray flux limits from the ROSAT All Sky Survey are taken from Table 8 of Plotkin 
et al. (2010a). 

The SEDs of the 34 objects are shown in Fig. lAll The majority of these sources have SEDs consistent with those of typical 
radio-quiet quasars in Richards et al. (2006), showing they are more likely to be WLQs in nature rather than BL Lac objects. 
However, there are also several sources with other kinds of SED profiles, which we discuss in more detail below. 

J0834+5112, J1556+3854, J1610+3039, and J1658+61 18 — These four sources have very red SEDs (strong in the infrared 
and weak in the optical/UV). They also have very red SDSS spectra. They are more likely to be absorbed quasars than bona fide 
WLQs. 

J 1421 +0522 — This source has an SED profile peaking in the near-infrared band, which is more similar to the SEDs of BL Lac 
objects. 

J1522+4137 — The SED of this source appears similar to that of J0945H-1009, which can be fit with a cold accretion disk 
model as in Laor & Davis (2009). It is probably a high-redshift, radio-quiet WLQ. 

J1633+4227 — The SED of this source peaks in the near-infrared band and drops rapidly in both blueward and redward 
directions. This SED profile is similar to that of a radio-loud BL Lac object with significant host-galaxy contamination (e.g., 
J0823H-1524 in Plotkin et al. 201 1). J1633H-4227 also has the same Ca II H/K break value (C = 0.33) as that of J0823H-1524. It is 
possible that J1633H-4227 is a radio-quiet BL Lac object with substantial host-galaxy contamination, although its radio-loudness 
needs to be constrained more tightly (currently (Xro> —0.19 estimated from its FIRST coverage). The near-infrared peak of the 
SED of J1633H-4227 could also perhaps be caused by extreme variability. 
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Fig. A1 . — Spectral energy distributions (SEDs) of additional radio-quiet WLQ candidates in Plotkin et al. (2010a) with "WISE coverage. All the lines and 
symbols follow the same definitions as in Fig.|9] The v-axis is in arbitrary units. All the SEDs are in the rest frame, except those for the two sources without 
redshift information (J1541+2631 and J1556+3854) which are plotted in the observed frame. 
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Fig. Al. — Continued. 
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TABLE 1 
X-RAY Observation Log 



Object Name (SDSS J) 




Aopt-x'' 
(arcsec) 


Detector 


Observation 
Date 


Observation 
ID 


Exp. Time 
(ks) 


Off-axis Angle 
(arcmin) 


References 


Chandra Cycle 12 Objects 


















noincr\ "7n i c^oCTi^ o 

Uoli5U. /y + 5zi5JU.5 


1.153 


0.8 


ACIS-S 


2010 Dec 28 


12710 


4.1 


0.3 


1 


uy45J3.y8 + luoysu. i 


1.671 




ACIS-S 


^A11 T 10 

201 1 Jan 12 


12706 


3.0 


0.3 


2 


110938.50 + 373611.7 


0.397 


0.3 


ACIS-S 


2011 Feb 27 


12711 


3.1 


0.3 


1 


125219.47 + 264053.9 


1.289 


0.3 


ACIS-S 


2011 Mar 12 


12709 


3.4 


0.3 


1 


153044.08 + 231013.4 


1.406 


0.4 


ACIS-S 


2011 Apr 15 


12707 


3.0 


0.3 


1 


161245.68 + 511816.9 


1.595 


0.4 


ACIS-S 


2011 Feb 01 


12708 


3.2 


0.3 


I 


Archival X-ray Data Objects 


















101353.46 + 492758.1 


1.640 




MOS' 


2004 Apr 23 


0206340201 


22.7 


6.4 


1 


113y00.55- 020140.0 


1.903 


0.2 


ACIS-S 


2004 Jul 21 


4871 


14.9 


0.6 


I 


160410.22 + 432614.6 


1.538 


0.2 


ACIS-I 


2006 Jun 25 


6933 


26.7 


3.7 


1 






0.2 


ACIS-I 


2006 Jun 23 


7343 


19.4 


3.7 




211552.88 + 000115.5 


2.500 




ACIS-S 


2008 Dec 24 


10388 


9.5 


0.3 


1,3,4 


232428.43+144324.3 


1.417 


0.7 


ACIS-S 


2009 May 31 


10386 


5.0 


0.3 


1,3,4 



References. — (1) Plotkin et al. 2010a; (2) Hryniewicz et al. 2010; (3) Collinge et al. 2005; (4) Plotkin et al. 2010b. 

Redshift for each source. See ^3.1 | for details about redshift measurements. 
^ Angular distance between the optical and X-ray positions; no entry indicates no X-ray detection. 
^ This object was observed by both the MOS and pn detectors. We list MOS detector parameters here. 



TABLE 2 

Quasar UV Emission-Line Measurements 



Object Name 


MJD 


C IV Blueshift 


REW 


REW 


REW 


REW 


REW 


(SDSS J) 






(CIV) 


(Si ivy 


(A 1900 A)'' 


(Fe III) 


(Mgll) 


Chandra Cycle 12 Objects 
















081250.79 + 522530.8 


53297 








3.8 + 2.1 


<4.5 


8.4 + 0.7 


094533.98+100950.1 


52757 


-7300+1700 


3.0+1.2 




4.9+1.5 


1.7+1.5 


17.1 + 0.8 


110938.50 + 373611.7 


53499 














125219.47 + 264053.9 


53823 








8.8+1.2 


2.9 + 0.9 


8.7 + 0.4 


153044.08 + 231013.4 


53878 








4.9+1.5 


2.6+1.2 


12.95 + 0.4 


161245.68 + 511816.9 


52051 


-4700+1300 


3.4+1.8 




5.1 + 1.5 


3.0+1.5 


9.5 + 0.6 


Archival X-ray Data Objects 
















101353.46 + 492758.1 


52076 








3.4+1.8 


<6.9 


6.2 + 0.8 


113900.55-020140.0 


52294 




<9.0 


<9.9 


< 10.8 


<9.0 


11.1 + 1.1 


113900.55-020140.0 (HET) 


55702 


-2950+1550 


3.2 + 2.7 




11.7+1.8 


5.5 + 1.5 




160410.22 + 432614.6 


52756 








< 1.9 


< 1.8 


5.8+1.0 


211552.88 + 000115.5 


52443 














232428.43+144324.3 


52258 






7.6 + 2.1 


<5.4 




8.6 + 0.8 


PHL 1811 " 




-1400 + 250 


4.7 + 0.9 


4.8 + 0.9 


8.3 + 0.6 


4.7 + 0.6 




vol composite '^•'^ 




-570+30 


30.0 + 0.3 


8.7 + 0.3 


21.7 + 0.2 


2.9 + 0.1 





Note. — The blueshift values are in units of km s ' . All REW values are in units of A. 
^ This line is a blend of Si IV and O IV]; we refer to it as Si IV simply for convenience. 

^ Mainly C III] A 1909, but also including [Neiii] A 1814, Sin A 1816, Aim A 1857, Sim] A 1892, and several Fe ill multiplets (see Table 2 of Vanden Berlc 
et al. 2001). 

^ These measurements are taken from Wu et al. (201 1). 
^ The composite spectrum from Vanden Berk et al. (2001). 



TABLES 

Continuum S/N Values for SDSS DR7 Quasars 



SNnao (1650-1750 kf 



SN3000 (2950-3050 A)"^ 



SNsiso (5100-5200 A)*^ 



Object Name (SDSS J) 


nii" 




MJD 


Plate 


Fiber 


S/N 


N - " 


1 e 
'^lo.obs 


2 t 

'^hi.obs 


S/N 




3 e 
'^lo.obs 




S/N 


N ■ " 

'•pix 


3 e 


J t 

'^hi.obs 


000006.53 + 003055.2 


20.09 


1.825 


52203 


685 


467 


4.20 


257 


4660 


4943 


2.22 


146 


8333 


8616 


0.00 











000008.13 + 001634.6 


19.49 


1.837 


52203 


685 


470 


4.83 


256 


4682 


4965 


3.87 


146 


8370 


8654 


0.00 











000009.26+151754.5 


19.15 


1.199 


52251 


751 


354 


2.30 


57 


3798 


3848 


10.97 


146 


6485 


6705 


0.00 











000009.38+135618.4 


18.30 


2.234 


52235 


750 


82 


14.22 


256 


5337 


5660 


0.00 











0.00 











000009.42- 102751.9 


18.77 


1.845 


52143 


650 


199 


9.64 


256 


4695 


4979 


5.61 


146 


8393 


8678 


0.00 












Note. — This table is available in its entirety in the online journal. Part of the table is shown here for guidance about its format. 
The apparent /-band magnitude using the BEST photometry of the SDSS DRV quasar catalog. 
The redshift value from the SDSS DR7 quasar catalog. 
Wavelength range in the rest frame. 
Number of pixels used in the calculation. 

The lower limit of the observed-frame wavelength range used in the calculation. 
The upper limit of the observed-frame wavelength range used in the calculation. 
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TABLE 4 

X-RAY Counts and Basic Spectral Properties 



Object Name (SDSS J) 


Full Band 
(0.5-8.0 keV)" 


Soft Band 
(0.5-2.0 keV)" 


Hard Band 
(2.0-8.0 keV)'' 


Band 
Ratio'' 




Chandra Cycle 12 Objects 

081250.79 + 522530.8 
094533.98+100950.1 


T ^+3.0 

< 3.3 


3 1+3-0 

^■'-1.7 

<3.2 


<3.3 
< 3.3 


< 1.06 


>0.72 



110938.50 + 373611.7 552.3+|^;^ 406.4+|i;| 145.8^[3;[ 0.36+°;g^ 1-68+°;^^ 

125219.47 + 264053.9 5.51^ ^■'^-Ij <7-l < 2-25 > -0.01 

153044.08 + 231013.4 I25.8+I7? 101.9+iJ[ 23.9+^° 0.23+j5J^^ 2.12+°^ 

161245.68 + 511816.9 39.3+^3' 30.5^1^ T-'^tH '^■^^-0.10 2.05+SJ;^§ 

Archival X-ray Data Objects 

101353.46 + 492758.1 < 38.6 <25.7 < 28.8 

113900.55-020140.0 5.3+3;^ 4.1+5"^ <5.1 < 1.25 > 0.55 

160410.22 + 432614.6 46.5lX^ 41.3+|;| 8.2l^;^ 0-201o;ot ^-^^-oH 

il.Stfl 26.3;fJ < 14.8 <0.57 > 1.25 

211552.88 + 000115.5 < 3.3' < 3.2 < 3.3 

232428.43+144324.3 6.4^3^ 41+3.2 2.I+2] 0.52+0;^; l-SS^ii' 

" Errors on the X-ray counts were calculated using Poisson statistics corresponding to the IcT significance level 
according to Tables 1 and 2 of Gehrels (1986). 

The band ratio is defined here as the number of hard-band counts divided by the number of soft-band counts. 
The en'ors on the band ratio coiTespond to the IfJ significance level and were calculated using equation (1.31) in 
§1.7.3 of Lyons (1991). The band ratios for all of the Chandra objects observed in the same cycle can be directly 
compared with one another. 

The effective power-law photon indices were calculated using the Chandra PIMMS tool (version 3.9i:). The 
effect of the quantum efficiency decay over time at low energies of the ACIS detector was corrected for Chandra 
observed objects. 



TABLE 5 

X-RAY, Optical, and Radio Properties 









Count 


logL 


log Ly 
(2500 A) 










Object Name (SDSS J) m," 


M, 




Rate'' 


^'0.5-2keV'^ /2keV'' (2-lOkeV) f^^^^^ 




Affox {of 






(1) (2) 


(3) 


(4) 


(5) 


(6) (7) (8) (9) 


(10) 


(11) 


(12) 


(13) 


(14) 



Chandra Cycle 12 Objects 



Archival X-ray 



081250.79 + 522530.8 


17.99 


-26.18 


4.33 


0.76^ 


-0.72 
-0.41 


0.34 


1.08 


43.48 


2.16 


30.88 




-2.03 


-0.42 (2.10) 


12.43 


> 


-0.093 


094533.98+100950.1 


17.44 


-27.54 


2.90 


< 1. 


,04 


<0.44 


< 1.77 


< 44.00 


3.42 


31.39 


< 


-2.03 


< -0.34 (2.32) 


> 7.69 


> 


-0.074 


110938.50 + 373611.7 


18.01 


-23.50 


1.57 


132.38 


.+6.84 
'-6.51 


54.61 


113.79 


44.55 


0.83 


29.42 




1.10 


0.31 (1.90) 


0.16 




-0.180 


125219.47 + 264053.9 


17.71 


-26.65 


0.75 


0.92^ 


-0.87 
-0.49 


0.37 


1.27 


43.64 


2.53 


31.04 




-2.03 


-0.39 (2.69) 


10.38 




-0.109 


153044.08 + 231013.4 


17.53 


-27.13 


4.23 


35.07; 


+3.75 
-3.40 


15.54 


55.76 


45.35 


3.32 


31.23 




-1.45 


0.22 (1.50) 


0.27 




-0.184 


161245.68 + 511816.9 


17.56 


-27.33 


1.66 


9.48^ 


-2.05 
-1.71 


3.92 


15.17 


44.89 


3.34 


31.34 




-1.67 


0.02 (0.11) 


0.89 


> 


-0.054 


ata Objects 


































101353.46 + 492758.1 


18.23 


-26.68 


0.79 


<0. 


,85 


<0.44 


< 1.72 


< 43.97 


1.59 


31.04 


< 


-1.91 


< -0.26 (1.81) 


>4.76 


> 


-0.115 


113900.55-020140.0 


18.88 


-26.43 


2.58 


0.27j 


-0.21 
-0.13 


0.12 


0.52 


43.55 


0.83 


30.88 




-2.00 


-0.38 (1.91) 


9.77 


> 


-0.172 


160410.22 + 432614.6'' 


17.84 


-26.94 


1.22 


1.56^ 


-0.21 
-0.19 


0.99 


3.74 


44.29 


2.21 


31.13 




-1.83 


-0.18 (1.22) 


2.94 


> 


-0.091 


211552.88 + 000115.5 


17.84 


-26.94 


6.17 


<0. 


,33 


<0.14 


< 0.73 


< 43.97 


2.39 


31.55 


< 


-2.12 


< -0.41 (2.78) 


> 11.70 




-0.138 


232428.43 + 144324.3 


19.22 


-26.81 


4.26 


0.82^ 


-0.64 
-0.39 


0.33 


1.18 


43.73 


1.21 


30.80 




-1.92 


-0.32 (1.60) 


6.82 


> 


-0.145 



^ The apparent /-band magnitude using the BEST photometi'y of the SDSS DRV quasar catalog. 
^ The count rate in the observed-frame soft X-ray band (0.5 — 2.0 keV) in units of 10~"^ s~^ . 

The Galactic absoiption-conected observed-frame flux between 0.5 — 2.0 keV in units of 10"'"^ ergs cm~" s~' . 
^ The flux density at rest-frame 2 keV, in units of 10~-^- ergs cm~- s~' Hz~^ 
^ The flux density at rest-frame 2500 A in units of 10~-^ ergs cm~" s~^ Hz~^ 

Aotox: the difference between the measured ctyx and the expected a^y., defined by the a^y. ~^^5QQ ^ relation in equation (3) of Just et al. (2007); the statistical significance of this difference, a, is measured in units of the RMS Ofox defined 
in Table 5 of Steffen et al. (2006). 

^ The factor of X-ray weakness compared to a typical radio-quiet quasar with similar optical/UV luminosity; see §3. 
^ The X-ray properties for J 1604-1-4326 reported here are for the average of its two Chandra observations. 
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TABLE 6 
Results of Peto-Prentice Tests 



Sample I (No. of sources) vs. 


Sample II (No. of sources) 


Statistic 


Null-hypothesis Probability 


Low-zRQWLQs(Il) 


RQ Sample B quasars (132) 


4.982 


6.29 X IQ-' 


High-z RQ WLQs (9) 


RQ Sample B quasars (132) 


2.693 


7.08 X 10-3 


Low-; + High-z RQ WLQs (20) 


RQ Sample B quasars (132) 


4.589 


4.45 X 10-'* 


Low-z RQ WLQs (11) 


High-z RQ WLQs (9) 


0.932 


0.351 



Note. — See Feigelson & Nelson (1985) for the detailed definition of the test statistic. The null-hypothesis 
probability was calculated from each test statistic using a Gaussian distribution, e.g., 1 — Pq = 6.29 x 10"', where 
Pq is the cumulative Gaussian probability at 4.982fT. 



TABLE 7 
X-RAY Spectral Analysis 





Power Law 
with Galactic Absorption 


Power Law 
with Galactic and Intrinsic Absorption 


Object Name 


r 




r Nh(\0^^ cm-2) 




Jl 109-1-3736 
J1530-H2310 
J1604+4326/J1612-H5118 


1 77+0.14 
'■"-0.14 

2 11 +0-37 

'^■'^"-0.30 


42.09/23 
14.18/9 
77.38/105" 


1.79+[;-^^ <0.14 
2.18+0;^] <2.67 
210^a33 <l-58 


42.08/22 
14.13/8 
77.38/105'' 



" The numbers here are C/n instead of j} l\, where C is the C-statistic, and n is the total number of spectral 
bins. 



